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Abstract 

(^ ' Using chiral perturbation theory we investigate the QCD shear viscosity (r/) to entropy density 

o' 

CN ' (s) ratio below the deconfinement temperature (~ 170 MeV) with zero baryon number density. 

It is found that rj/s of QCD is monotonicahy decreasing in temperature (T) and reaches 0.6 

lO ' with estimated ~ 50% uncertainty at T = 120 MeV. A naive extrapolation of the leading order 

(N 

result shows that rj/s reaches the 1/4tt minimum bound proposed by Kovtun, Son, and Starinets 



> 

oo 
en 



using string theory methods at T ~ 200 MeV. This suggests a phase transition or cross over 
might occur at T < 200 MeV in order for the bound to remain valid. Also, it is natural for r]/s 



(^ \ to stay close to the minimum bound around the phase transition temperature as was recently 

^^ ■ found in heavy ion collisions. 
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I. INTRODUCTION 

Shear viscosity rj characterizes how strongly particles interact and move collectively 
in a many body system. In general, strongly interacting systems have smaller rj than 
the weakly interacting ones. This is because r] is proportional to eTmft, where e is the 
energy density and Tmft is the mean free time, which is inversely proportional to particle 
scattering cross section. Recently a universal minimum bound for the ratio of rj to entropy 
density s was proposed by Kovtun, Son, and Starinets |il|. The bound, 

^ > - , (1) 

s Att 

is found to be saturated for a large class of strongly interacting quantum field theories 
whose dual descriptions in string theory involve black holes in anti-de Sitter space [2, y. 



Recentfy^/s close to the minimum bound were found in relativistic heavy ion collisions 
(RHIC) p, llllsl. This discovery came as a surprise. Traditionally, quark gluon plasma 
(QGP) — the phase of QCD above the deconfinement temperature Tc{^ 170 MeV at zero 
baryon density J2|) — was thought to be weakly interacting. Partly because lattice QCD 
simulations of the QGP equation of state above 2Tc were not inconsistent with that of 
an ideal gas of massless particles, e = 3p, where e is the the energy density and p is 
the pressure of the system J9| . However, recent analyses of the elliptic flow generated by 
non-central collisions in RHIC JT], |8[ and lattice simulations of a gluon plasma [lO| yielded 
rj/s close to the the minimum bound at just above X- This suggests QGP is strongly 
interacting at this temperature.* (However, see Ref. [l^ for a different interpretation.) 

Given this situation, one naturally wonders if rj/s of QCD was already close to the 
minimum bound at just above Tc, what would happen if we keep reducing the temperature 
such that the coupling constant of QCD gets even stronger? Will the t]/s minimum bound 
hold up below T^? If the bound does hold up, what is the mechanism? Is the change of 
degrees of freedom through a phase transition or cross over sufficient to save the bound? 
If the bound does not hold up, what is the implication to string theory? 

To explore these issues, we use chiral perturbation theory (xPT) and the linearized 
Boltzmann equation to perform a model independent calculation to the rj/s of QCD in 

* See als o llll Il2l Il3| . For discussions of the possible microscopic stucture of such a state, see [ij, |la, 
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the confinement phase. Earher attempts to compute meson matter viscosity using the 
Boltzmann equation and phenomenological phase shifts in the context of RHIC hydrody- 



namical evolution after freeze out can be found in Refs. J24 , 



m 



26| . In the deconfinement 



phase, state of the art perturbative QCD calculations of r] can be found in Refs. |22l, |2 

II. LINEARIZED BOLTZMANN EQUATION FOR LOW ENERGY QCD 

In the hadronic phase of QCD with zero baryon-number density, the dominant degrees 
of freedom are the lightest hadrons — the pions. The pion mass 171,^ = 139 MeV is much 
lighter than the mass of the next lightest hadron — the kaon whose mass is 495 MeV. Given 
that Tc is only ~ 170 MeV, it is sufficient to just consider the pions in the calculation of 
thermodynamical quantities and transport coefficients for T <^ T^.. 

The interaction between pions can be described by chiral perturbation theory (xPT) in 



Ltion tneory 
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xPT is a 



a systematic expansion in energy and quark {u and d quark) masses 
low energy effective field theory of QCD. It describes pions as Nambu-Goldstone bosons 
of the spontaneously broken chiral symmetry. At T ^ T^, the temperature dependence 
in vrvr scattering can be calculated systematically. At T = T^, however, the theory breaks 
down due to the restoration of chiral symmetry. ^" 

The shear viscosity rj of the pion gas can be calculated either using the Boltzmann 
equation or the Kubo formula. Since the Boltzmann equation requires semi-classical 
descriptions of particles with definite position, energy and momentum except during brief 
collisions, the mean free path is required to be much greater than the range of interaction. 
Thus the Boltzmann equation is usually limited to low temperature systems. The Kubo 
formula does not have this restriction. In this approach rj can be calculated through the 
linearized response function 

r} = -\j dt' f dt /"dx3([r^(0),r^(x,t)]) (2) 

with T*-' the spacial part of the off-diagonal energy momentum tensor. One might think 
a perturbative calculation of the above two point function will give the answer for rj. But 
this can not be true if r] a Tmft-, as mentioned above, for T„ift — > 00 in the free case. Indeed, 



^ The QCD chiral restoration temperature and the deconfinement temperature happen to be close to 
each other at zero baryon density. We do not distinguish the two in this paper. 



the Kubo formula involves an infinite number of diagrams at the leading order (LO) 
However, in a weak coupling 0^ theory, it is proven that the summation of LO diagrams 



is equivalent to solving the linearized Bo^ 
particle masses and scattering amplitudes 



tzmann equation with temperature dependent 
23l | . This proof extended the applicable range of 
the Boltzmann equation to higher temperature but is restricted to weak coupling theories. 
In the case we are interested (QCD with T < 140 MeV), the pion mean free path is always 
greater than the range of interaction (~ 1 fm) by a factor of 10^. Thus, even though the 



3, 



the temperature is still low 



coupling in xPT is too strong to use the result of Ref. 
enough that the use of the Boltzmann equation is justified. 

The Boltzmann equation describes the evolution of the isospin averaged pion distri- 
bution function / = /(x, p,t) = fp{x) (a function of space, time and momentum) as 

^dMx) = ^ / dTu;3p {/l/2(l + /3)(1 + fp) - (1 + /l)(l + /2)/3/p} , (3) 

^p ^ J 123 

where Ep = -^/p^ + m^and g^^ = 3 is the degeneracy factor for three pions , 

and where T is the scattering amplitude for particles with momenta 1,2^ 3, p. In xPT, 
the LO isospin averaged vrvr scattering amplitude in terms of Mandelstam variables (s, t, 
and u) is 

\r\' = \ Y. (2/+l)|r(^f = ^{21mt + 9.^-24M> + 3(t-nf} . (5) 

^ 7=0,1,2 ^^-^ 

The temperature dependence in pion mass and pion scattering amplitudes can be treated 
as higher order corrections. 

In local thermal equilibrium, the distribution function fp (x) = (^e^'^^)^i^^^)p^ — ij 
with i3{x) the inverse temperature and V^{x) the four velocity at the space-time point x. 
A small deviation of fp from local equilibrium is parametrized as 

fp{x) = fl,'\x)[l-{l + fl,'\x)}xp{x)] , (6) 

while the energy momentum tensor is 



(27r)3 Ep 



T^u{x) = g^ / 'L'^ fpix) . (7) 



We will choose the V(x) = frame for the point x. This implies d^V^ = after tak- 
ing a derivative on V^{x)V^{x) = 1. Furthermore, the conservation law at equilibrium 
9^T^'^|j^p=o = allows us to replace dtP{x) and dtY{x) by terms proportional to V ■ V(a;) 
and V/3(a;). Thus, to the first order in a derivative expansion, Xpi^) can be parametrized 

as 

X,ix) = Pix)Aip)V-Vix)+Pix)B{p) (p,,p, - ^S,,^ (^YiYM±1j}M _ l^^^.v . V(x)^ 

(8) 
where i and j are spacial indexes. "^ A and B are functions of x and p. But we have 
suppressed the x dependence. 

Substituting (jH)) into the Boltzmann equation, one obtains a linearized equation for B 

IpiPj - i^^ijP^ j = ^V^ / c?ri2;3p(l + ni)(l + n2)n3(l + np)"^ 

X [Bij (p) + Bij ih) - Bi,{k2) - % (fci)] = gJi, [B] , (9) 
Bij{p) = B{p) ipiPj - -6,j 

where we have dropped the factor (ViVj{x) + VjV^(x) — trace) contracting both sides of 
the equation and write /^ (x) at this point as rij = (e^'^' — l) . There is another integral 
equation associated with V ■ V(x) which is related to the bulk viscosity ( that will not 
be discussed in this paper. The V ■ /? and (9iV terms in p^d^fp will cancel each other by 
the energy momentum conservation in equilibrium mentioned above. 

In equilibrium the energy momentum tensor depends on pressure V{x) and energy 
density e(x) as T^J{x) = {V{x) + e{x)} V^{x)V„{x) — V{x)6^„. A small deviation away 
from equilibrium gives additional contribution to T^i, whose spacial components define 
the shear and bulk viscosity 

ST.j = -2V ^YjyM±lj}M. _ ^5^^.v . V(x)) + CS.,V ■ V(x) . (10) 

After putting everything together we obtain 
g^p f d^p 



V 



-Up (1 + Up) Bij{p) (piPj - g^iiPM 



10 J (27r)3Ep 
^ / I^^^^P (1 + np) B,,{p)F,, [B] ^ gl{B\F[B]) . (11) 



■f A non-derivative term is not allowed since fp should be reduced to fp when /3 and V^ become 
independent of x. There is no term with single spacial derivative on f3{x) either. The only possible 
term (V • V) /3(x) vanishes in the V(a;) = frame. 



Here we see immediately that if the scattering cross section is scaled by a factor a, 

c^ri2;3p -> a ((iri2;3p) , (12) 

then Eqs. (jHl) and pij) imply the following scaling 

Bij{p) -> a'^Bij{p) , 

V -^ "~^^ , (13) 

with 7] proportional to the inverse of scattering cross-section. This non-perturbative result 

is a general feature for the linearized Boltzmann equation with two-body elastic scattering. 

To find a solution for B{p), one can just solve Eq. Q. But here we follow the approach 



outlined in Ref. 



25 



26| to assume that B{p) is a smooth function which can be expanded 
using a specific set of orthogonal polynomials 

oo 

B{p) = g-'\p\yY,brB'^^\4p))^ (14) 

r=0 

where B^^\z) is a polynomial up to z^' and br is its coefficient. The overall factor |p|^ will 
be chosen by trial and error to get the fastest convergence. The orthogonality condition 

/ ^^E;""' ^^ + ""'^ \p\yB(^\z)B(^\z) oc 6r,s (15) 

can be used to construct the B'^^'\z) polynomials up to normalization constants. For 
simplicity, we will choose 

B^'\z) = l. 

With this expansion, the consistency condition for B{p) in Eq. ljlll) yields 

r/ = ^6,L« = 5^6,(i?W|F[i?W])6, , (16) 

r r,s 

where 

^^'^ = TeI ^^W""' ^^ + ""'^ \p\'B^'\p) « V • (17) 

Since br is a function of m^r, f-^ and T, the fo^'s in Eq. (ll6p are in general independent 
functions, such that L*^*"^ = '^^{B^^^F[B^^'>])bs [one can show that this solution of bg gives 
a unique solution of rj] , or equivalently 

V^W = 5^(5M|F[5(^)])6, . (18) 



This will allow us to solve for the bg and obtain the shear viscosity 

rj = boL^^\ (19) 

In the next section, we will show that this expansion converges well, such that one 
does not need to keep many terms on the right hand side of Eq. pSj) . If only the s = 
term was kept, then 

n ~ W^ . (20) 

(5(o)|F[5(o)]) 

The calculation of the entropy density s is more straightforward since s, unlike r], does 
not diverge in a free theory. In xPT, the interaction contributions are all higher order in 
our LO calculation. Thus we just compute the s for a free pion gas: 

s = -9.P ^-^ , (21) 



where the partition function Z for free pions is 

\ogZ If d=^p 



up to temperature independent terms. 

III. NUMERICAL RESULTS 

In this section we present the results for t] and r]/ s of QCD up to T = 120 MeV at zero 
baryon number density. In Fig. 1 the LO xPT result of r] using the linearized Boltzmann 
equation is shown. The lines with circles, squares and triangles correspond to keeping 
the first one, two and three polynomials on the right hand side of Eq. ()18p . respectively. 
We have used y = and z{p) = |p| to construct the polynomials. The figure shows the 
expansion converges rapidly. As a test of the calculation, we also reproduce the shear 
viscosity result of Ref. [23j for 0^ theory by setting the scattering amplitude T = A to 
be a constant. In 0^ theory, rj is monotonically increasing in T. If T ^ m^, rj oc T^/X^ 
with T^ given by dimensional analysis and A^^ by the scaling of coupling shown in Eqs. 
(I12|l and ()13|1 . In xPT, however, rj is not monotonic in T. At T ^ 171^^, the scattering 
amplitude is close to a constant, thus xPT behaves like a 0^ theory. But at T ^ ?7i^, T 
(X T'^/ f^ and rj oc f^/T. At what temperature the transition from rj cc T^ to rj cc 1/T 
takes place depends on the detail of dynamics. In xPT, this temperature is around 20 
MeV. 
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FIG. 1: (Color online) Shear viscosity as a function of temperature in LO xPT. 1st, 2nd and 
3rd are the results of keeping the first one, two and three polynomials on the right hand side of 
Ea. (|18j) . respectively. 

The radius of convergence in momentum for xPT is typically 47r/7r ~ 1 GeV. To 
translate this radius of convergence into temperature, we compute the averaged center 
of mass momentum (|p|) = ^J {B\p'^\F[B]) / {B\F[B]) . We found that for T = 120 and 
140 MeV, (IpI) ^460 and 530 MeV < Anf^. However, xPT is expected to break down 
at the chiral restoration temperature (~ 170 MeV). Thus our LO xPT result can only 
be trusted up to T ~ 120 MeV. At the next-to- leading order (NLO), it is known that 
the isoscalar vrvr scattering length will be increased by ~ 40% [22]. This will increase the 
cross section by ~ 100% and reduce ?7 by ~ 50% near threshold. This is an unusually 
large NLO correction since a typical NLO correction at threshold is < 20%. The large 
chiral corrections does not persist at the higher order. At the next-to- next-to-leading 
order (NNLO), the correction is ~ 10% ^]. Thus, to compute r] to 10% accuracy, a NLO 
xPT calculation is needed. 

The LO xPT result for r]/s is shown in Fig. 2 (line with rectangles). The error is 
estimated to be ~ 50% up to 120 MeV. ri/s is monotonically decreasing and reaches 0.6 
at T = 120 MeV. This is similar to the behavior in the m^r = case (shown as the line 
with rectangles) where rj/s oc f^/T"^ with s oc T^ from dimensional analysis and /tt = 87 
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FIG. 2: (Color online) Shear viscosity to entropy density ratios as functions of temperature. 
Line with circles (rectangles) is the LO xPT result with rriTr = 139 (0) MeV and fj^ = 93(87) 
MeV. Line with triangles is the result using vrvr phase shifts (PS). Dashed line is the conjectured 
minimum bound I/Att ~ 0.08. 

MeV in the chiral limit 

For comparison, we also show the result using phenomenological vrvr phase shifts |29| for 
1] but free pions for s. (Our result for rj is in good agreement with that of 26| for T between 
60 and 120 MeV. For an earlier calculation using the Chapman-Enskog approximation, 
see Ref. |30[.) This amounts to take into account part of the NLO tttt scattering effects 
but ignore its temperature dependence and the interaction in s. Since not all the NLO 
effects are accounted for, this t]/s is not necessarily more accurate than the one using LO 
xPT. The comparison, however, gives us some feeling of the size of error for the LO result 
we present here. Thus, an error of ~ 100% at T = 120 MeV for the LO result might be 
more realistic. 

Naive extrapolations of the three rj/s curves show that the 1/Att = 0.08 minimum 
bound conjectured from string theory might never be reached as in phase shift result (the 
first scenario), or more interestingly, be reached at T ~ 200 MeV, as in the LO xPT 
result (the second scenario). In both scenarios, we see no sign of violation of the universal 
minimum bound for t]/s below Tc. But to really make sure the bound is valid from 



120 MeV to Tc^ a lattice computation as was performed to gluon plasma above Tc [lO[ is 
needed. In the second scenario, assuming the bound is valid for QCD, then either a phase 
transition or cross over should occur before the minimum bound is reached at T ~ 200 
MeV. Also, in this scenario, it seems natural for r]/ s to stay close to the minimum bound 
around T^ as was recently found in heavy ion collisions. 

In the second scenario, one might argue that the existence of phase transition is already 
known, otherwise we will not have spontaneous symmetry breaking and the corresponding 
Nambu-Golstone boson theory at low temperature in the first place. Indeed, it is true 
in the case of QCD. However, if the rj/ s bound is really set by Nature, then a phase 
transition is inevitable in the vicinity of the temperature where the bound is reached. For 
a spontaneous symmetry breaking theory, the general feature of rj/ s we see here seems 
generic. At very high T, collective motion is weak, thus rj/ s gets smaller at lower T. At 
very low T in the symmetry breaking phase, the Nambu-Goldstone bosons are weakly 
interaction at low temperature, thus r]/ s gets smaller at higher T. A phase transition 
should occur before the extrapolated rj/ s curve coming from high T reaches the bound at 
Ti. Similarly, a phase transition should occur before the extrapolated rj/ s curve coming 
from low T reaches the bound at T2. Thus the range of phase transition is Ti < Tc < T2. 
However, it is also possible that the first scenario takes place and r]/ s bounces back to 
higher values without a phase transition. In this case, it is less clear what makes r]/ s 
non-monotonic and it certainly deserves further study. 

It is interesting to note that the degeneracy factor (7^ drops out of r] while the entropy 
s is proportional to g.,, as in Eqs. ffTHjl and (pTjl . respectively. This suggests the 77/5 bound 
might be violated if a system has a large particle degeneracy factor.^ For QCD, large g.,^ 
can be obtained by having a large number of quark flavors Nf with g^^ ~ iV?. However, 
the existence of confinement demands that the number of colors N^ should be of order Nf 
to have a negative QCD beta function. After using /^^ oc y/Wc) the combined N^. and Nf 

scaling of tj/ s is 

r? f^ N'^ 

which is of order one. Thus QCD with large Nc and Nf can still be consistent with the 
Tjl s bound below T^. 



We thank Thomas Cohen for pointing this out to us. This possibihty was also mentioned in Ref. [l|. 



IV. CONCLUSION 

We have explored whether the conjectured t]/ s minimuni bound will hold up below the 
QCD deconfinement temperature. Using chiral perturbation theory and the linearized 
Boltzmann equation we have computed the QCD r]/ s ratio at zero baryon number density 
and for T < 120 MeV. It is found that t]/s is monotonic decreasing in T and it reaches 
0.6 with estimated 50% uncertainty at T = 120 MeV. Naive extrapolations have shown 
that rj/s met the l/Air minimum bound conjectured from string theory at T ~ 200 
MeV as in the LO xPT case, or rj/s stayed above the bound as in the phenomenological 
phase shift case. In the former case, in order for the t]/s lower bound to remain valid at 
higher temperature, a phase transition or cross over should occur at T < 200 MeV before 
the bound is reached. We argued that this might be a general feature for spontaneous 
symmetry breaking theories that the extrapolation of the low(high) temperature rj/s curve 
sets a upper(lower) bound on Tg. Our result also suggests that it is natural for rj/s to 
stay close to the lower bound around the phase transition temperature as was recently 
found in heavy ion collisions. 

As this paper was being finished, reference J22] appeared. In that paper, some of the 
relations between Tc and the rj/s bound are also discussed. 
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